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Abstract. The XMM-Newton spectra of the isolated neutron star RX J0720.4— 3125 show deviations from a Planckian en- 
ergy distribution below 400 eV, similar to the spectra of RBS1223, another long-period X-ray pulsar, as reported recently by 
iHaberl et al. (.2003) . For a Gaussian-shaped absorption line we derive an energy of 271 eV and an equivalent width of —40 
eV from the phase-averaged spectra of RX J0720.4— 3125. We investigate the spectral variations seen in hardness ratios as 
function of pulse phase and find that they are best described by changes in the depth of the absorption line. The line equivalent 
width changes between —31 eV around intensity maximum of the pulse and —58 eV at the declining part of the pulse. Small 
variations (<20 eV) of the line energy with pulse phase may still be caused by statistical fluctuations. On the other hand, the 
black-body temperature varies significantly by 2.5 eV (statistical 90% errors typically 0.7 eV) reaching the highest value at 
pulse maximum. One possible interpretation for the absorption line is cyclotron resonance scattering of protons in a magnetic 
field with B : 
pulsar 



: 5 X 10 G. This field strength is compatible with estimates inferred from recent spin down measurements of the 
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1. Introduction 

The soft X-ray source RX J0720.4— 3 125 was dis covered 
in the ROSAT all-sky survey data by iHaberl et alJ ^997). 
The detection of 8.39 s sinusoidal pulsations in the X- 
ray flux, its black-body-like spectrum with little interstel- 
lar a bsorption and the ident i fication with a faint blue optical 
star (Motch & Haberl 1998; Kulkar ni & van Kerkwiiklll998t 



iKaplan et al..2003;,Motch et al..2003.) classified this object as a 
nearby isolated neutron star. It is the second brightest of a group 
of seven, maybe eight, radio quiet (recently stringent limits 
were obtained for RX J0720.4-3125 and RX J0806.4-4123, 
|johnstonll2003l) and X- ray dim isola ted neutron stars (XDINs; 
for a recent review see Haberl 2003'). 

The available low-resolution X-ray spectra obtained by 
ROSAT from XDINs were all consistent with Planckian en- 
ergy distributions of temperatures kT in the range 40 - 100 
eV and little attenuation by interstellar absorption. More re- 
cent observations of RX J0720.4— 3 125 and the brightest object 
RX J 1856.4—3754 were performed using the low energy trans- 
mission grating spectrometer (LETGS) aboard Chandra and 
the reflection grating spectrometers (RGSs) of XMM-Newton. 
Also at high spectral resolution both sources show no signifi- 
cant narrow features in their spectra which can best be modeled 
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I2OOII) or a modified Planckian spectrum (IBurwitz et alJl2003t) . 

In particular the high signal to noise LETGS spectrum from 
a 500 ks Chandra observation of RX J1856.4-3754 did not 
reveal any significant deviatio n from a smooth continuum 
spectrum ( Bur witz et alJ l2003l) . The statistical quality and 
energy band coverage of the LETGS and RGS spectra of 
RX J0720.4— 3125 are, however, insufficient to detect broad 
and shallow features in the spectrum. 

Cyclotron resonance absorption features in the 0.1-1.0 keV 
band can be expected in spectra from magnetized neutron stars 
with field strengths in the range of lO^O-lO" G or 2x10^3- 
2x10^'* G if caused by electrons or protons, respectively. The 
X-ray pulsar RBS1223 was the first XDIN for which a sig- 
nificant narrow-band deviation (at ~300 eV, equivalent width 
-150 eV, width ^100 eV) fr om a Planckian spectrum was dis- 
covered taaberl et al."2003'. hereafter H03). This absorption 
feature was discussed in terms of proton cyclotron resonance 
scattering in a magnetic field of (2-6) x 10^'^ G. 

H03 also suggested that changes in the soft part of the X- 
ray spectru m of the pulsar RXJ07 20.4— 3125 with pulse phase 
reported by Cropper et al. (200 1', hereafter COl), are caused 
by variable cyclotron absorption. In this paper we present the 
results from a spectral analysis of XMM-Newton data on the 
isolated neutron star RX J0720.4— 3125. The observed spectra 
show a broad absorption feature which varies with pulse phase. 
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Table 1. XMM-Newton EPIC observations of 
RXJ0720.4-3125. 



Time(UT) XMM Read-out Count rate 
Start End Detector Mode'^' [s^^]^^^ 


Exp. 

[ks] 


2000 May 13, Satellite Revolution 078, Obs. ID 0124100101 


02:43 19:59 MOSl/2 FF/SW 1.63/1.87 
02:30 20:00 pn FF 7.59 
01:42 15:36 RGSl/2 SQ 0.26/0.23 


62.0 
62.5 
50.0 


2000 Nov. 21-22, Satellite Revolution 175, 0132520301 


21:18 02:21 MOSl/2 LW 1.31/1.40 
19:21 02:36 pn FF 6.51 
18:30 03:06 RGSl/2 SQ 0.24/0.23 


18.2 
26.1 

30.9 


2002 Nov. 6-7, Satellite Revolution 533, 0156960201 


17:52 02:14 MOSl/2 FF 

18:14 02:15 pn FF 7.66 

17:51 02:16 RGSl/2 SQ 0.23/0.22 


30.0 
28.4 
30.2 


2002 Nov. 8-9, Satellite Revolution 534, 0156960401 


19:25 04:18 MOSl/2 FF 

19:47 04:18 pn FF 7.64 

19:24 04:19 RGSl/2 SQ 0.23/0.22 


31.8 
30.2 
32.0 



FF: Full Frame (time resolution of 0.073 s and 2.6 s for EPIC- 
pn and -MOS, respectively; SW: Small Window (MOS time resolu- 
tion 0.3 s); LW: Large Window (MOS time resolution 0.9 s); SQ: 
Spectro-l-Q 

Mean count rates are for the energy bands 0.13-1.5 keV (EPIC) 
and 0.32-0.90 keV (RGS); when not given the data are not used in the 
spectral analysis. 

The thin filter was used in all EPIC cameras in revolutions 78, 533 and 
534 and the medium filter in revolution 175. 



We discuss this in the framework of proton cyclotron resonance 
or atomic line transitions in a strong magnetic field. 

2. XMM-Newton observations 

The X-ra y pulsar RX J0720.4 -3 125 was observed with XMM- 
Newton (Jansen et al. 2001) on five occasions. Here we uti- 
lize the data collected with the European Photon Imaging 
Cameras (EPICs) which provide medium energy-resolution X- 
ray spectra. The c ameras are based on MOS (EPIC- MOS 1 
and -M OS2. lTurner et al. 2001) andpn (EPIC-pn. Struder et alJ 
I2001I) CCD technology and are mounted in the focal plane of 
the three X-ray multi-mirror systems ( lAschenbac h et al. 2000). 
Simultaneous high spectral resolution data were obtained with 
the Reflection Grating Spectrometers dden Herder et al ■I2OOII) . 
The details of the XMM-Newton observations are summarized 
in Table [2 We do not use the last observation because EPIC- 
pn was operated in SW read-out mode (less well calibrated at 
low energies) and with the thick filter During the first obser- 
vation the EPIC instruments suffered from strong background 
flares which were screened out leaving only 36 ks (pn) and 39 
ks (MOS) which were used for the spectral analysis. From the 
other three observations all the data was used. 



In processing and preparation of the observational data and 
model fitting of the spectra we followed H03 (using an extrac- 
tion radius of 40" for the EPIC data). We used the same version 
of the XMM-Newton analysis software SAS and calibration 
files. 

3. Pulse-phase averaged spectra 

3.1. Combined XMM-Newton spectra 

H03 reported the analysis of the XMM-Newton spectra of 
RBS1223. They could not find satisfactory fits using a sim- 
ple black-body model, non-magnetic neutron star atmosphere 
models with different element abundances or two black-body 
components with different temperatures. An acceptable fit was 
found with a broad absorption line included in the black- 
body model. In the following analysis of the spectra of 
RX J0720.4— 3125 we find very similar residuals when com- 
paring the measured spectra to the simple absorbed black-body 
model. We therefore also include a Gaussian-shaped absorp- 
tion line in the modeling of the spectra of RX J0720.4— 3125. 
To account for cross-calibration problems between the different 
instruments we allow the black-body temperature to be fit indi- 
vidually to each spectrum but fit the absorption column density 
as single free parameter to all sixteen XMM-Newton spectra 
simultaneously. Similarly, all three parameters of the absorp- 
tion line (energy, width and normalization) were fit to all spec- 
tra simultaneously. The results for both models are reported in 
Table 12] The line at 271 eV is relatively broad with a width 
of a = 64 eV. The line equivalent width (EW) is about —40 
eV', less prominent as measured in the spectra of RBS1223 
with EW = -150 eV (H03). The XMM-Newton spectra of 
RX J0720.4— 3125 together with the best fit models are shown 
in Fig.[n 

The deviation from a Planckian shape is highly significant 
as can be seen from the improvement in the values (Table|2} 
when an absorption line is included ( Ax^=683 for 3 additional 
parameters in the fits to the spectra of the EPIC-pn instrument 
which is most sensitive at energies below 300 eV). In Fig. |2l 
the incident photon spectrum is shown (for clarity only the 
summed EPIC-pn spectrum from the three thin filter observa- 
tions is plotted). This presentation better visualizes the absorp- 
tion feature as compared to Fig. [0 where the detector response 
at energies below 500 eV (with relatively bad energy resolu- 
tion and strong energy redistribution into lower detector chan- 
nels due to electron losses in the surface layers of the detector) 
smears out the intrinsically broad feature. 

3.2. A problem in the spectrai calibration? 

IPaerels etaP fEoOll presented a black-body fit to the EPIC-pn 
spectrum of RX J0720.4— 3125 obtained from the observation 
in May 2000. Residuals around 300 eV were attributed to un- 
certainties in the filter transmission data. However, the model- 
ing was based on early spectral response calibration data and 

' All values for EW are obtained with the spectral fitting package 
XSPEC version 11.2.0. 
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Fig. 1. Simultaneous fits using models A (left) and B (right) to the XMM-Newton spectra of RX J0720.4-3125. For model 
definition see Table |2] For each model the best fit (histogram) to the spectra (crosses) is plotted in panels (a). Panels (b) - (d) 
show the residuals for EPIC-pn, -MOS and RGS spectra, respectively. For model B panel (e) illustrates the best fit model with 
the absorption line removed. The three EPIC-pn spectra obtained with thin filter were combined for clarity in the plots, as well as 
all the eight RGS spectra. The MOS data below 300 eV were not used for the spectral fits. The residuals increasing with energy 
above 800 eV in the EPIC spectra are probably caused by pile-up (see Sect. 13. 3> . 



RXJ0720.4-3125 was the only XDIN observed with XMM- 
Newton at that time. In addition to RX J0720.4-3125 and 
RX J1856.4-3754 four XDINS are now observed with XMM- 
Newton, partly several times and different instrumental setup 
and we have access to all the data. The spectral modeling of the 
spectra is in progress and revealed significant deviations from 
a Planckian shape for RBS1223 (H03) RX J0720.4-3 125 (as 
presented here) and RX J1605.3H-3249 jvan KerkwiilJ2004 . 

The best quality X-ray spectrum of any XDIN was obtained 
from RXJ1856.4-3754 with the Chandra LETGS. It is well 
described by a Planckian modeP and can be used as a calibra- 
tion reference for the soft X-ray band. Black-body fits to the 
EPIC-pn spectrum of RX J1856.4— 3754, using the same in- 
strument calibration as for RX J0720.4— 3125 are fully accept- 



^ As discussed inl Burwitz et alj i2003h the broad-band spectrum of 
RX J1856.4— 3754 is even better modeled by a modified black-body 
spectrum of the form x Bi, (B,j is the Planck function). We note 
that the spectral deviations of RX J0720.4— 3125 discussed in this pa- 
per are found in a relatively narrow energy band in contrast to the 
smooth wide-band modification caused by the power-law component 
which mainly changes the black-body parameters. 



able, with deviations less than 5% at all energies below 1 keV-'. 
For clarity and better comparison with RX J0720.4— 3125 we 
show the EPIC-pn spectrum of RX J1856.4-3754 with best 
black-body fit in Fig. |3l The residuals in the black-body fit 
to the spectra of RX J0720.4-3125 (Fig. □ clearly show a 
strong systematic behavior, which is not seen in the case of 
RX J 1856.4— 3754. Hence, we conclude that the deviations 
from a Planckian shape revealed in the EPIC-pn spectra of 
RX J0720.4— 3125 are not caused by calibration uncertainties. 
In addition, the strength of the observed feature varies with 
pulse phase (see Sect. 4). 

However, as can be seen from Table|2] there are still cross- 
calibration problems between the different instruments which 
cause systematic uncertainties larger than the statistical errors. 
The EPIC-MOS spectra are barely consistent to each other at 
low energies and there seems to be a systematic difference be- 
tween the residuals of the model fit for thin and medium filter 
spectra (medium filter spectra better agree with those of EPIC- 

' See the presentation by M. Kirsch at 
http://xmm.vilspa.esa.es/extemal/xmm_user rsupport/usersgroup7l 
2003033 1/presentations.shtml 
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Table 2. Spectral analysis of pulse-phase averaged X-ray spectra of RX J0720.4— 3 125. 





Model A: phabs*bbody^^' 






Model B: phabs*(bbody+gaussian)''^^ 




Obs-Inst. 


kT 


Nh 


x'/dof 


kT 


Nh 


Eline 


Cline 




X^/dof 






[eV] 


[10^°cm-2] 


per inst. 


[eV] 


[10^° cm"^] 


[eV] 


[eV] 


[eV] 


per inst. 


[erg cm~^ s~^] 


078-pn 


03.2 


1 10 

L.JQ 




04 3-1-0 3 


79±0 1 3 


T71 1 1 A 


o4±/ 


— j9.0 




1.19x10"^^ 


1 /^-pn 


84.6 


_(4) 








— 


— 


—39.4 




1 99 V 1 n~^^ 

1 .ZZ. /s. J.U 




88.6 


— 




07 <-\-0 A 




— 


— 


—40.6 




1 . 1" A iAj 




88.3 


— 


1679/639 


07 O-i-O zL 




— 


— 


—40.5 


996/636 


1 IQvIO"^''^ 
1 . 1 7 X lU 


VJ / o-iVHJo 1 


84.4 


— 




8A 1 -1-0 A 




— 


— 


—39.6 




1 91 V 1 (\~^^ 




83.8 


— 








— 


— 


—39.4 




1 . ID X lU 


078 Mn*\9 

u / iviwoz. 


81.6 


_- 








— 


— 


—38.8 




1 AO V 1 {\~^^ 

1 X lU 


1 / J-WIkJoZ. 


82.7 


— 


408/167 


89 A-l-0 8 




— 


— 


—39.1 


468/167 


1 98 V 1 n~i^ 

i .Zo X lU 


078 Rri*v 1 


86.7 






8'^ O-l- 1 1 








—39.8 




QQ w 1 n~^i 


1 7^^ Rn*vi 


89.6 






87 74-1 
/ . / Ill l.J 




z 


Z 


-40.7 




QzL V 1 n ~ ^ 1 
u.y'+x iU 


533-RGSl 


95.0 


= 




93.1ibl.6 


— 


= 


= 


-42.7 




0.86x10-^^ 


534-RGSl 


94.1 






92.2±1.6 


= 






—42.4 




0.88x10"^^ 


078 Rri*v9 


86.7 






8S ^-1- 1 1 








-39.9 




QA V 1 n~i^ 


175-RGS2 


88.4 






87.0±1.4 








-40.5 




0.92x10"" 


533-RGS2 


92.8 






91.5±1.4 








-42.1 




0.85x10"" 


534-RGS2 


94.4 




1453/1266 


92.9±1.4 








-42.6 


1455/1266 


0.85x10"" 


total 






3540/2072 












2919/2069 




LETGS 


82.5 


1.50 


157/181 


78.8±2.0 


1.56±0.20 


27 If 


64f 




157/181 


0.99x10"" 



All errors are given for a 90% confidence level. Parameters fixed in the fit are marked by "f behind the value. 
Standard XSPEC names for the model components. 

For MOS and RGS spectra the line EW is calculated from the extrapolated model because of the restricted Une coverage. 
Observed flux 0.1-2.4 keV. 
(4) (ienotes fit parameter is forced to be the same for all spectra. 




0.2 0.5 1 

Energy (keV) 



Fig. 2. Incident ('unfolded') photon spectrum for the black- 
body model with absorption line. The crosses represent the 
combined EPIC-pn spectrum from the three thin-filter obser- 
vations. The best fit black-body model with absorption line and 
with the line normalization set to zero are plotted as solid and 
dashed histograms, respectively. 




0.1 0.2 0.5 1 

Channel Energy (keV) 



Fig. 3. Black-body model fit to the EPIC-pn spectrum of 
RXJ1856.4— 3754. The lower two panels show the residuals 
from the best fit in units of a and as ratio data/model. 



pn down to 200 eV). We excluded therefore the MOS data be- 
low 300 eV in the spectral analysis. The EPIC-pn spectra are 
in large agreement with each other, although those obtained 
from the last two observations yield significantly higher black- 



body temperatures. These temperature differences are caused 
by apparent gain variations by up to 3% at energies below 1 
keV which were measured in emission line spectra of super- 
nova remnants (routine calibration targets to monitor the en- 
ergy scale) and the instrumental Oxygen K-edge energy deter- 
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mined from strong sources with flat continuum spectra . The 
origin of this effect is not well understood and can presently not 
be corrected for. Gain variations shift the spectral part above 
500 eV which results in the temperature change, while the spec- 
trum below 500 eV, which determines the Nh (and where the 
deviations from a black-body model are seen), is practically 
unaffected (allowing Nh to vary in the EPIC-spectra yields val- 
ues consistent within the errors). To account for the gain varia- 
tions we therefore allowed kT to vary between the spectra. The 
somewhat higher flux^ derived from the revolution 175 spec- 
trum may be caused by reduced pile-up losses (see Sect. 13. 3> . 

3.3. Pile-up in tlie EPIC spectra 

Bright sources observed with X-ray CCD detectors can be sub- 
ject to pile-up effects. If the read-out is not sufficiently fast, 
two or more photons may deposit charge in one or few neigh- 
boring CCD pixels during the exposure of a single read-out 
frame. The energy of the detected event will then be the sum of 
the energies of the original photons. Therefore pile-up hardens 
a continuum spectrum. RX J0720.4— 3 125 was observed by the 
EPlC-pn in full frame mode which provides pile-up free spec- 
tra for count rates below ~6 cts s^^. As can be seen from the 
count rates listed in Table[2the EPlC-pn spectra are somewhat 
affected by pile-up. This probably explains the residuals in the 
model fits at energies above ~800 eV (Fig.0. Also the EPlC- 
pn spectrum obtained with medium filter (which reduces the 
total count rate) shows the lowest pn-derived black-body tem- 
perature and the highest flux (lowest flux losses). 

Pile-up occurs predominantly in the center of the point 
spread function (PSF) where the probability is highest that two 
photons are recorded in the same pixel within a single read- 
out frame. The SAS task 'epatplot' can be used to visualize the 
pile-up effects on the spectrum. The main effect on very soft 
spectra like those of XDlNs is a spectral hardening which leads 
to a somewhat higher black-body temperature and to some 
flux-loss when using single-pixel events only. In the case of 
RX J0720.4-3 125 we excluded the inner area of the PSF with 
20" in radius (15" for the medium filter observation with lower 
count rate) to obtain practically pile-up free spectra. The dis- 
advantage of this procedure is in reduction of the number of 
extracted counts by a factor of 7.3 (4.0 in the medium filter 
case). The statistics becomes almost insufficient to investigate 
the deviations from a black-body spectrum and therefore, we 
fixed the line energy and width at the values derived in Table|2] 
Table|3lsummarizes the black-body model parameters infeiTed 
from the pile-up free EPlC-pn spectra. The derived tempera- 
tures are typically 3 eV lower than those obtained from the 
spectra not corrected for pile-up effects. The column density 
derived for both cases is consistent within the errors. There is 
again a significant difference in temperature between the first 
and last two spectra. However, at the cuiTent status of the cali- 
bration it is not clear if this is a long-term systematic effect. 



See the presentation by F. Haberl at 
^tp://epic3. xra.le.ac.uk/pub/cal-pv/meetings/tuebingen-20 03-02l 
^ Errors on the pn fluxes are 7x 10"^'* erg cm~^ s~^ 



007B thin 
0175 medium 

0533 thin ^i^^ 

0534 thm ^fe 









0.1 0.3 0.5 1 

Energy (keV) 



Fig. 4. EPlC-pn spectra extracted from a ring-shaped area (in- 
ner and outer radii are 20"and40" for the thin filter, and 15"and 
40" for the medium filter observations), excluding the inner 
part of the PSF in order to avoid pile-up effects. 



Table 3. Fits to pUe-up free EPlC-pn spectra. 



Model B: phabs*(bbody+gaussian) 








Obs-Inst. 


kT 


Nh 


Elinc 


Clinc 


x'/dof 




[eV] 


[10^°cm-^] 


[eV] 


[eV] 


per inst. 


078-pn 


81.6±0.9 


0.90±0.13 


27 If 


64f 




175-pn 


81.6±0.9 










533-pn 


85.2±1.0 










534-pn 


84.7±1.0 








477/439 



3.4. Comparison with other instruments 

Chandra observed RX J0720.4-3125 in February 2000 with 
the LETGS for a total of -38 ks (observation IDs 368, 369 
and 745). The LETGS spectrum was found to be consistent 
with a black-body model with kT ::^80 eV and absorp tion col- 
umn density -1.7x10^" cm"^ jPavlov et al]l2002l) . We re- 
analyzed these data and obtained kT = 82.5 ±4.0 eV and Nh = 
(1.5±0.5)x 10^" cm^^ for a single black-body model. To ver- 
ify that the LETGS spectrum is also consistent with the pres- 
ence of an absorption feature we added a Gaussian-shaped line, 
with parameters as derived from the XMM-Newton spectra, to 
the black-body model. The fit results for both models, with 
and without line, are listed in Table |2 The statistical quality 
of the LETGS spectrum is clearly insufficient to discriminate 
between the models (Fig.|5jl. The column density inferred from 
the combined fit to the XMM-Newton spectra is significantly 
lower than the LETGS value. Allowing individual values of 
Nh for the three XMM-Newton instruments (all other param- 
eters in model B are treated as described in Table |2} results 
in Nh/10^° cm-2 = L23±0.07, 0.66±0.30 and 0.73±0.30 for 
EPlC-pn, -MOS and RGS, respectively, revealing still cross- 
calibration uncertainties. The EPlC-pn value is in best agree- 
ment with the LETGS value which is probably because the 
LETGS spectrum of RX J1856.4-3754 was used for calibra- 
tion of the low-energy spectral response of EPlC-pn. 
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Fig. 5. Fits using a black-body model with (bottom pan- 
els) and without absorption line to the LETGS spectra of 
RX J0720.4— 3125. The low statistical quality of the spectrum 
does not allow to discriminate between the models. 

4. Pulse-phase variations 

Because of the highest sensitivity, better time resolution 
(Tabled and full line coverage we investigate only EPIC-pn 
spectra in our pulse-phase spectroscopic study. In the 0.12- 
1.2 keV energy band the pulse profile of RX J0720.4-3125 
is nearly sinusoidal. We fitted a sine wave to the folded light 
curves obtained from the EPIC-pn data in this band to deter- 
mine the pulsed fraction. An example is shown in Fig.|6l The 
pulsed fraction (after background correction) is around 11%. 
Pulse periods and pulsed fractions inferred from the four EPIC- 
pn data sets are listed in Table 0] The presented pulse peri- 
ods are used here to extract the phase-resolved s pectra and we 
refer f or a more detailed temporal analysis to Cropper et al. 

The smaller pulsed fraction found in the revolution 
175 data collected with medium filter (which more strongly 
absorbs low-energy photons), suggests an energy dependence 
of the pulse profile. To investigate this further we determined 
the pulsed fraction as a function of energy as plotted in Fig.E] 
We combined the data from the three thin filter observations 
assuming a period change of 5 x 10^^"' s s^^ cZane et al..2002) 
when folding the data. There is an indication that the pulsed 
fraction is larger below '^400 eV compared to that at higher 
energies. We note here that due to the strong re-distribution 



Table 4. X-ray pulsations derived from the EPIC-pn data of 
RXJ0720.4-3125. 



Rev. 


Pulse 


Pulsed'^' 




Period 


Fraction 




[s] 


[%] 


078 


8.3911165 


10.5±0.3 


175 


8.3912190 


10.1±0.4 


533 


8.3911276 


11.1±0.4 


534 


8.3912461 


11.1±0.4 



0.12-1.2 keV 




0.5 1 l.D 2 

Phase (B. 3912461 s) 



Fig. 6. Pulse profile of RX J0720.4-3125 obtained from the 
observation in revolution 534 together with the best fit sine 
wave. 

at low energies the data points below 400 eV are not statisti- 
cally independent (as discussed in Sect. 3.1 and demonstrated 
in panel "e" of Fig. [Q. Also, this approach is not sensitive to 
changes in the folded light curves related to the shape of the 
pulse profiles. In the next step we therefore regard only two 
energy bands below and above 400 eV. In Fig.|8lthe two light 
curves are plotted together with their ratio (hardness ratio). A 
dependence of the hardness ratio on pulse phase was first re- 
ported by COl. At energies above 400 eV the pulse peak be- 
comes broader and more flat-topped. This results in a hardness 
ratio minimum near the intensity maximum and two (asymmet- 
ric) hardness ratio maxima at the flanks of the intensity peak. 
A test for constancy of the hardness ratio yields a reduced 
value xt = 4.6 (for 29 degrees of freedom). The corresponding 
probability that the hardness ratio is constant is extremely low: 
2.5x10-15. 

We divided the pulse into four phase intervals of equal 
length in order to investigate the X-ray spectrum as function 
of pulse phase. The intervals cover the intensity minimum, the 
rise, the intensity maximum and the decline of the pulse. The 
latter phase interval is centered at the hardness ratio maximum 
derived from the EPIC-pn data. Compared to the approach 
above, spectral modeling has the advantage that the spectral 
response (energy resolution, re-distribution) is taken into ac- 
count. However, for statistical reasons only a small number 
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Fig. 7. Pulsed fraction as function of energy obtained from the 
combined EPIC-pn thin filter data of RX J0720.4-3125 from 
revolutions 78, 533 and 534. For each energy bin a light curve 
with 10 phase bins was folded and the pulsed fraction deter- 
mined by (Nt-NDc)/Nt with Nt beeing the total number of 
counts and Ndc the number of counts at the DC level of the 
light curve. The energy bins were chosen to obtain at least 5000 
events for each phase bin. 

of phase bins can be investigated. The sixteen phase-selected 
spectra (from four observational data sets) were fit again simul- 
taneously with the model including an absorption line. Because 
of the reduced statistics the line width was fixed at 64 eV as ob- 
tained from the phase-averaged spectra. 

To verify if the hardness ratio variations can be caused by 
changes in the absorption line parameters we fit the column 
density simultaneously as a single free parameter to all six- 
teen spectra. The energy of the line (kept equal for the spec- 
tra from the same phase interval) and the black-body temper- 
ature (for all spectra) were allowed to be fit individually. The 
resulting fit parameters are summarized in Table |5] The spec- 
tra from "soft" and "hard" phases (as defined in Table |5} are 
presented in Fig.|5l No significant shift of the line energy be- 
tween the two phase intervals is seen. The depth of the line 
varies strongly, which results in an EW during the "soft" phase 
(around intensity maximum) about a factor of two smaller than 
that during the "hard" phase (at the end of the declining part of 
the pulse). At the same time, there is a 2-3 eV increase in the 
black-body temperature found during the intensity maximum. 
At this phase the softening due to reduced line absorption dom- 
inates the spectral hardening caused by increased temperature 
and yields the hardness ratio minimum. 

To see if the variations in the line depth are really required 
to account for the spectral changes, we forced the line depth 
to be a single fit parameter, kept the same for all four phase 
intervals. This results in a fit with x^/dof = 2394/1970 in com- 
parison to 2267/1967 with free variation in line depth (Table|5}. 
Comparison of the two models using an F-test yields a proba- 
bility of 4.3 X lO^^'^ that the line depth is not variable. Allowing 
variations in column density Nh instead of variable absorption 




Phase (8.39112 s) 



Fig. 8. Pulse profile of RX J0720.4-3125 obtained from the 
combined EPIC-pn thin filter data as used in Fig. in two en- 
ergy bands, together with the hardness ratio. The line in the 
bottom panel marks the best-fit with a constant hardness ratio 
which yields a x^=4.6 for 29 dof. 

line (depth and energy are constant with phase) yields x^/dof = 
2336/1970 with reduced Nh during "soft" phase as was already 
reported by COl. Also for this case the F-test probability of 
9.7 X 10^^"^ strongly favors the model with variable line depth. 
Our new results show that the model with phase-dependent ab- 
sorption line yields the best description of the phase-resolved 
spectra of RX J0720.4-3125. 

5. Discussion 

From a thorough spectral analysis of four XMM-Newton ob- 
servations of the isolated neutron star RX J0720.4— 3125 we 
find deviations in the spectra from a Planckian shape, similar 
to those reported by H03 for RBS1223. A model consisting of 
a black-body continuum attenuated by photo-electric absorp- 
tion and an absorption line yields acceptable fits to the XMM- 
Newton spectra and is also consistent with the Chandra LETGS 
spectrum currently available from RX J0720.4— 3 125. From the 
pulse-phase averaged spectra we derive a line energy of 271 
eV, (T = 64 eV and EW = —40 eV, for a Gaussian-shaped line. 
These absolute values may still be subject to systematic uncer- 
tainties larger than the statistical errors due to cross-calibration 
problems of the different instruments used for our analysis, but 
relative changes of model parameters, e.g. with pulse phase, 
could be determined to high accuracy. 

Hardness ratio variations with the pulse phase from 
RX J0720.4— 3125 were akeady reported by COl, indicating 
spectral variations. When formally fit with an absorbed black- 
body model the spectrum near intensity maximum with lower 
hardness ratio yields a reduced Nh compared to phases around 
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Table 5. Pulse-phase resolved spectra. 



Model B: phabs*(bbody+gaussian); x^/dof = 2267/1967 



Rev. kT 


Nh 


Elinc 


Clinc 


EW 


[eV] 


[10 cm J 


[eV] 


[eV] 


[eV] 


pulse minimum 










078 83.0±0.6 


0.88±0.10 


285±6 


64f 


-46.1±3.9 


175 82.4±0.7 


= lW 


= 1 


= 1 


-46.0±3.9 


533 86.5±0.7 




= 1 


= 1 


-47.1±3.9 


534 86.2±0.7 




= 1 


= 1 


-47.0±3.9 


pulse rise 










078 84.2±0.6 




288±8 


= 1 


-31.6±4.2 


175 83.3±0.7 




=5 


= 1 


-31.5±4.2 


533 88.0±0.7 




=5 


= 1 


-32.4±4.2 


534 87.0±0.7 




=5 


= 1 


-32.2±4.2 


pulse maximum, "soft" phase 








078 84.2±0.6 




271 ±9 


= 1 


-3 1.1 ±4.4 


175 83.2±0.7 




=9 


= 1 


-30.9±4.4 


533 87.2±0.7 




=9 


= 1 


-31.9±4.4 


534 87.2±0.7 




=9 


= 1 


-31.9±4.4 


pulse decline, "hard" phase 








078 82.4±0.5 




279±5 


= 1 


-57.3±3.5 


175 81.4±0.6 




=13 


= 1 


-57.0±3.5 


533 84.9±0.6 




=13 


= 1 


-58.3±3.5 


534 84.9±0.6 




=13 


= 1 


-58.3±3.5 



"=n" denotes fit parameter is linked with parameter in line n 



hardness ratio maximum. COl suggested that this could be ex- 
plained by energy-dependent beaming effects. Alternatively, 
they discuss variable absorption by matter captured by and 
rotating with the magnetosphere of the neutron star. Taking 
into account the absorption line revealed in the phase-averaged 
spectra we find that the modulation of the hardness ratio with 
pulse phase can be better explained by variations in the absorp- 
tion line depth which lead to changes in the line EW by about 
a factor of two. During pulse rise and at intensity maximum, 
where we see a hardness ratio minimum, the absorption line is 
weakest. The black-body temperature varies by 2-3 e V between 
"hard" and "soft" phase spectra, with the highest temperature 
during intensity maximum. The resulting hardness ratio mini- 
mum shows that the decrease of the absorption line depth is the 
dominating effect near pulse maximum. 

H03 suggested cyclotron resonance absorption as a likely 
origin for the absorption feature seen in the X-ray spectra of 
RBS1223. Since charged particles as electrons, protons or ions 
can be the origin of this absorption one cannot uniquely de- 
rive the magnetic field strength from the cyclotron line energy 
without additional independent information. First estimates for 
the magnetic field strength of RX J0720.4— 3 125 were infeiTed 
from measurements of the s pin down ra t e unde r the assumption 
of magnetic dipole braking lzane et al give a magnetic 

field B^IO^'^ G, which would exclude electrons as the origin 
of the cyclotron line, which is expected at energies around 100 
keV for such a strong field. In the case of protons a cyclotron 
line energy of 271 eV would result in a magnetic field strength 
of 5.3 X lO^'^ G for a neutron star with canonical mass and ra- 



spectrum of the X-ray pulsar RXJ0720.4-3 125 
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Fig. 9. Phase-resolved EPIC-pn spectra from phases of high 
and low hardness ratios (see Fig. fig-hardness). The upper pair 
of spectra shows the combined data from the three thin filter 
observations, the lower pair corresponds to the medium filter 
observation. The upper and lower spectra in each pair are ex- 
tracted from the phases of low ("soft" phase) and high ("hard" 
phase) hardness ratio. 

dius (1.4 Mq, 10 km). Highly ionized atoms of heavy elements 
have a mass to charge ratio of ^^2 with respect to protons and 
would lead to B a factor of ~2 higher than that derived for pro- 
tons. Different ionization states would result in a series of lines 
with energies differing by only a few percent. Smearing out 
of these lines due to the B field variations in a rotating dipole 
field may cause a single broad absorption f eature, broader than 
that ex pected from the B field smearing teavlinet alJll995t 

IZane et al 2001 ) of a single line and hardly resolvable. In 

a most recent re-analysis of the availa ble archival and new 
XMM-Newton data of RX J0720.4 - 3 1 25 lCropper etaP t004 
constrain B to (2.8 - 4.2) x lO^'^ G from the pulse period history 
by interpreting the spin-down rate in terms of magnetic dipole 
radiation losses. This result is most compatible with the cy- 
clotron line scenario involving protons. 

An alternative possibility for the origin of the absorp- 
tion line is atomic bound-bound transitions. In strong mag- 
netic fields atomic orbitals are distorted into a cylindrical 
shape and the electron energy levels are similar to Landau 
states, with binding energies of atoms significantly increased. 
E.g. for hydrogen in a magnetic field of the order of 10^^ 
G the strongest atomic transition is expected at energy 
E/eV » 75(l-H0^131n( Bi.s))-H63Bi.-,, with B13 = B/IO^^ G 
JZavlin & Pavlovl2002l) . For the line energy found in the spec- 
tra of RXJ0720.4-3125 this would require B ~ 3x10^3 q 
Atomic line transitions are expected to be less prominent at 
higher temperatures because of a higher ionization degree 



Haberl et al.: A phase-dependent absorption line in the spectrum of the X-ray pulsar RXJ0720.4— 3125 



9 



o 

CD 



o 

LO 



o 



62 82.5 63 83.5 84 84.5 

kT (eV) 

Fig. 10. Black-body temperature vs. absorption line equivalent 
width derived from the EPIC-pn spectra from the four pulse 
phase intervals described in the text and listed in Table|5l(only 
the values derived from the first observation are plotted). 

llZavhn & Pavlovl i2002). Qualitatively, this is in line with the 
observed anti-correlation between temperature and equivalent 
width (see Fig.llO>. However, the observed effect is much too 
strong: a temperature change of 2 eV at a temperature of 83 
eV would not cause a change of the number of ionization states 
by a factor of two. In RBS 1223 the effect is even stronger 
as it shows an EW of -150 eV at similar temperature. At the 
same time in this source the pulsed fraction is much larger than 
in RX J0720.4— 3125 suggesting that the change of angles be- 
tween the magnetic axis and the line of sight is much larger. 
We conclude that the anti-correlation between temperature and 
equivalent width does not reflect a direct physical connection, 
but is due to the fact that both correlating quantities are depend- 
ing on the pulse phase. 

Therefore, we prefer the interpretation of the observed line 
feature in terms of the proton cyclotron resonance. The ob- 
served dependence of temperature and equivalent width on 
pulse phase may be due to the change in the viewing geometry 
of the inclined magnetic rotator. However, a thorough quantita- 
tive modeling of cyclotron and atomic bound-bound absorption 
effects in strong magnetic fields (a few lO^'^ G) of isolated neu- 
tron stars and their observable rotational phase dependences is 
required to understand the behavior in detail. 
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